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Abstract. This paper presents a new methodology to count the number of numerical semigroups of given 
• genus or Frobenius number. We apply generating function tools to the bounded polyhedron that classifies 

C ~ ) i the semigroups with given genus (or Frobenius number) and multiplicity. First, we give theoretical results 

about the polynomial-time complexity of counting the number of these semigroups. We also illustrate the 
methodology analyzing the cases of multiplicity 3 and 4 where some formulas for the number of numerical 
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1. Introduction 

A numerical semigroup is a subset S of N that is closed under addition, G S and generates Z as a 
group. This last condition is equivalent to gcd(5) = 1. 

For a given numerical semigroup S, the set G(S) — N\S, known as the set of gaps of S, has finitely 
many elements. Furthermore, S has a unique minimal system of generators {ni < ... < n p }. The 
element ri\ is the least positive integer belonging to S and it is denoted by m(S'), the multiplicity of S 
(1x1(5) = min(5 \ {0})) and the cardinality of G(S) is known as the genus of S, g(5). The largest integer 
■ not in S is known as the Frobenius number of S and it is denoted by F(5). The interested reader is 

00 • referred to [5D] for further details on numerical semigroups. 

£SJ . Given neS \ {0}, the Apery set of S with respect to n is the set Ap(5*, n) = {s G S : s — n $ S} and 

it can be easily shown that if for every i £ {0, . . . , n — 1} we take w(i) the least element in S congruent 
with i modulo n (denoted w(i) = i (mod n)), then Ap(5, n) = {0 = w(0),w(l), . . . ,w(n — 1)}. The set 
Ap(5, n) completely determines S, since S = (Ap(5, n) U {n}) (where (A) denotes the monoid generated 
by A). Moreover, the set Ap(S,n) contains, in general, more information than an arbitrary system of 
generators of S. For instance, Selmer in [21] gives the formulae, g(S) = (^w>eAp(s m(S)) w ) — 
^ ■ and F(S) = max(Ap(5, n)) — n. Moreover, for all s G S there exist unique t G N and w G Ap(5, n) such 

that s = tn + w. The smallest Apery set is Ap(S f , m(S)). 

Rosalcs ct al. |16j and previously Kunz in [14] gave a one-to-one correspondence between the set of 
numerical semigroups with multiplicity m and the set of integer points inside a rational polyhedron in 
R m_1 . This correspondence is based on the Apery set description of these numerical semigroups. These 
polyhedra are, in general, not bounded, but performing adequeate cuts to them, like fixing the gender or 
the Frobenius number, we obtain a polytope whose integer points are identified with each of the numerical 
semigroups with multiplity m and fixed gender or Frobenius number, respectively. Then, the problem 
of counting numerical semigroups is equivalent to the problem of counting the number of integer points 
inside a polytope. However, it is well-known that the problem of detecting a lattice point in polyhedra 
is NP-hard 10 . In this paper, we apply some results on short generating functions to provide some new 
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complexity results about the task of counting the number of numerical semigroups, and also some explicit 
descriptions for these numbers for multiplicities 3 and 4. 

A few papers have appeared recently analyzing the number of numerical semigroups with given genus 
or Frobenius number [51 [TH] . Here we give new results on the complexity of counting the number of 
numerical semigroups when the genus or the Frobenius number are fixed using short generating functions 
. A constructive approach for counting numerical semigroups and some conjectures about these numbers 
have been recently presented in [9] and [5j. In those approaches, the computation of the set of numerical 
semigroups with genus g requires to compute previously those with genus g — 1, proceeding recursively. 
Hence, giving exponential-time algorithms in g and consuming a considerable amount of CPU memory. 
However, although the proposed methodology here has nice theoretical properties, from a computational 
viewpoint the methods given in [9] or [5] are capable to compute all the numerical semigroups for genus up 
to 50, while we were only able to compute the number of numerical semigroups with genus up to 15 (see 
Table [TJ . One of the advantages of our methodology is that for given multiplicity, the polytope obtained 
fixing the genus (respectively the Frobenius number) can be seen as a parametric polytope, and a quasi- 
polynomial description of the number of numerical semigroups with any genus (respectively Frobenius 
number) can be computed. 

The paper is organized as follows. Section [2] recalls the main notions and results on rational generating 
functions. Section [3] is devoted to prove the complexity results about counting the number of numerical 
semigroups, maximal embedding dimension numerical semigroups with given genus, and a table with the 
results for genus up to 15. We also provide in this section a complexity result about computing the 
number of numerical semigroups with given Frobenius number and multiplicity. In Section [4] we give 
explicit formulas for the number of numerical semigroups and maximal embedding dimension numerical 
semigroups with multiplicities 3 and 4, fixing genus and Frobenius number with multiplicities three and 
four. 



2. Short generating functions 

Short generating functions were used by Barvinok [3], initially as a tool for counting the number of 
integer points inside convex polytopes. This tool is based in the geometrical papers by Brion [7] , Khovanskii 
and Puhlikov [11], and Lawrence [15] . 

Let P = {x € R" : A x < b} be a rational polyhedron in R™. The main idea is to encode the integer 
points inside a rational polytope as a "long" sum of monomials: 

f(P;z)= ]T z a , 
aePnZ" 

where z a = z" 1 ■ ■ ■ z" n . 

Barvinok's aimed goal was representing that formal sum of monomials in the multivariate polynomial 
ring Z[zi, . . . ,z n ], as a "short" (polynomially indexed, for fixed n) sum of rational functions. Actually, 
in [5], Barvinok provides a polynomial-time algorithm when the dimension, n, is fixed, to compute those 
functions. 
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Theorem 2.1 (Theorem 5.4 in [5]). Assume n, the dimension, is fixed. Given a rational polyhedron 
P C M™ , the generating function f(P; z) can be computed in polynomial time in the form 

z Ui 



/(P;z) = 5>- 



i=i 

where I is a polynomial- size indexing set, and where e £ {1, —1} and Ui,Vij G Z" for all i and j. 

As a corollary of this result, Barvinok gave an algorithm for counting the number of integer points in 
P. It is clear from the original expression of f(P; z) that this number is f(P; 1), but 1 = (1, . . . , 1) is a 
pole for the rational function. Hence, the number of integer points in the polyhedron is Km f(S; z). This 
limit can be computed using residue calculation tools from elementary complex analysis. 

A new algorithm for counting the integer points inside convex polytopcs using short generating functions 
was recently developed by Verdoolaege and Woods [25] and implemented in the software barvinok }22j . 
This software also allows to count integer points inside parametric polytopes |23j 

The above approach, apart from counting lattice points, has been applied to solve Operations Research 
problems. Actually, integer programs [3 [24], multiobjective integer programs [4], bilevel integer pro- 
grams [13] or integer programming games [12], among many others problems arc studied using Barvinok's 
rational functions. Furthermore, enumerating the complete set of solutions of a diophantinc system of in- 
equalities (in the bounded case) may be seen as multicriteria integer problem, and then, a polynomial-delay 
algorithm (for fixed dimension) have been developed in [4] for this task. 

3. Using Barvinok's algorithms to count numerical semigroups with given genus 

In this section we present some results concerning the task of counting the number of semigroups with 
fixed multiplicity and genus or Frobenius number. These results are based on the transformation of that 
problem to the problem of counting the number of integer points inside certain convex polytope. This is 
when Barvinok's rational function theory comes into scene. 

For the sake of completeness, we include the following technical result that we use in Theorem 13. II 

T)( X ) 

Lemma 3.1. Let g(x) = , where p(x) and q(x) are polynomials in the indeterminate x, and such that 

q{x) 

1 is a root of q(x) with multiplicity r. If lim g(x) exists, then p k '(l) = for k = 1, . . . , r. 

x^l 

Proof. Note that since 1 is a root of q(x) with multiplicity r, wc can write q(x) = (1 — x) r h(x) with h(x) 
a polynomial with h{\) ^ 0. 



Assume that p '(1) ^ for k < r. Clearly, q '(I) = 0, and then, lim g(x) = lim 



l x^i q(x) 

p k \x) p k Hl) 

lim , = = oo, a contradiction with the existence of limit of g(x) at x = 1. □ 

x— >i q k ) (x) 

Theorem 3.1. For fixed multiplicity m, counting the numerical semigroups for any genus g is doable in 
polynomial time. 

Proof. Let S(m) be the set of all numerical semigroups with multiplicity m 6 N\ {0, 1}. In [16] , Rosales et 
al. proved that there is a one-to-one correspondence between this set and the set of non-negative integer 
solutions of a system of linear diophantine inequalities. This identification was previously used by Kunz 
in HI. 
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Let S be in S(m) with Ap(5, m) = {0 = w(0), w(l), . . . , w(m - 1)}. For all i G {1, . . . , m - 1}, let 
fcj G N be such that w(i) = kim + i. Then (fci, . . . , k m -i) is a non-negative solution of the system 

Xi ^ 1 for all i G {1, . . . , m — 1}, 

(1) Xi + Xj — Xi+j for all 1 s; i j ^ m — 1, i + j ^ m — 1, 

+ aij — Xi+j- m —1 for all 1 < i < j < m — 1, i +j > m, 

Xi G Z for all i G { 1 , . . . , m — 1 } . 

Denote by T(m) the set of non-negative solutions of |T]). Then, T(m) and <S(m) are one-to-one identified 
(Theorem 11 in [16]). 

Selmer in [3T] proved that if Ap(5, m) = {0, w(l) = k\m + 1, . . . , w(n — 1) = fc ro -i + m — 1}, g(5) = 



5(SJ (E M eAp(S,m(S)) w ) - ^T 1 - and thcn thcn 5 = 51 fci 



Thus, the following system describes completely the family 5 m>g , the set of numerical semigroups with 
multiplicity m and genus g. 

Xi ^ 1 for all i G {1, . . . , m — 1}, 

Xi + Xj ~ Xi + j for all 1 ^ i ^ j ^ m — 1 , i + j ^ m — 1 , 

(2) + Xj — Xi+j- rn ^ —1 for all 1 ^ i ^ j ^ m — 1, i + j > m, 

m — 1 

i=l 

Xi G Z for all i G {1, . . . , m — 1}, 

Furthermore, @ defines a polytope (bounded polyhedron) in R m_1 , that we will call P m ,g- 

Each element in P m ,g H Z Tn_1 corresponds with a numerical semigroup with multiplicity m and genus 

g. We denote n m , 3 = #(P m , 9 D Z™- 1 ). 

By Theorem 12. 1[ the short generating function encoding the integer points inside P m ,g is computable 

in polynomial time for fixed m (that in this case is redundant since m is clearly fixed). Let f(P m ,g] z) = 
e,; — be that generating function. Choose c G Z m such that cuy ^ for all i,j (see [3] for 

further details about the polynomial-time complexity of this choice). Do the changes Xi = t Ci for each 

i = 1, . . . , m. Then, after those changes, let f(t) := f(P m ,g',t) = £ i~ • 

tEl 11(1 - t c ^) 

P(i) 

Clearly, we can write f(t) = , where P(t) and Q(t) are polynomial with <3(1) ^ 0. By 

(1 £)' Qit) 

using L'Hopital Rule (it can be done by Lemma IXTj) sequentially to compute the limit: 

lim f(P m . g ; z) = lim f(P m , g ; t) = lim J^^W) = ^T) = 

These operations are clearly polynomially bounded, and then, the overall procedure runs in polynomial 
time. □ 



For a numerical semigroup with fixed genus g, its multiplicity is at most g + 1 (this case is achived at the 
numerical semigroup (g+ 1, . . . , 2g + 1). Then, the number of semigroups with fixed genus g (independently 
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9+1 

of the multiplicity) is given by the finite sum n g = rim,g- By applying Theorem 13.11 a + 1 times we 

m=2 

obtain the following result. 

Corollary 3.1. Let g be a fixed genus, counting the number of numerical semigroups of genus g is doable 
in polynomial time. 

We have implemented the steps in the above proof to obtain the number of numerical semigroups of 
genus up to 15. Table [T] shows the number of numerical semigroups from genus 2 to 15 obtained using 
the above methodology We run the software barvinok [22| for counting n m>g , for each pair (m,g) with 
2<TO<<7 + lina PC with an Intel Pentium 4 processor at 2.66GHz and 1 GB of RAM. It was able 
to compute the number of numerical semigroups for genus up to 15. Although this methodology can be 
applied to any genus, the software barvinok fails to compute the short generating function representation 
of the polytope in ^ for g > 15. barvinok is designed to be applicable to general polytopes and therefore 
does not exploit the special structure of the polytope @. Up to date, there are only few implementations 
for computing generating functions of rational polyhedra. barvinok seems to be the most recent software 
and incorporates a large battery of options. More effective implementations oriented to the particular 
polytope in @ and further research in the computation of generating functions of polyhedra would permit 
go beyond the case g = 16. It is clear that, although theoretically our methodology is better than the 
complete enumeration proposed in [5l [6] , these other methods are faster than the one proposed here (the 
package numericalsgps [5] takes 0.08 seconds to compute the whole set of numerical semigroups up to 
genus 15 while our algorithm takes 38.85 seconds to count the set of numerical semigroups with multiplicity 
9 and genus 12). It seems to be due to the fact that the current implementation in barvinok does not 
take into account the special structure of the polytopes describing numerical semigroups. Next we apply 
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Table 1. Number of numerical semigroups with given genus g and multiplicity m. 



a similar methodology but fixing the Frobenius number instead of the genus, that is, our goal here is to 
count the number of numerical semigroups with fixed multiplicity and Frobenius number. 

Let S be a numerical semigroup with multiplicity m and Frobenius number F. The key here is that the 
Frobenius number of S can be computed by using the Apery set as F(S) = max Ap(S', m) — m (see j21j). 
Then, adding to the system of diophantinc inequalities |T]) those related to the Frobenius number: 

mxi + i < F + m, 

and the one that fixes the element in the Apery set that reaches the maximum: mxk* + k* = F + m (where 
k* = F (mod m)), we have a system of inequalities describing a polytope that characterizes completely 
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the set of numerical semigroups of multiplicity m and Frobcnius number F. 

Xi 1 for alH G { 1 , . . . , m — 1 } , 

Xi + Xj — Xi+j for all 1 ^ i ^ j ^ m — 1 , i + j ^ m — 1 , 

(3) Xi + xj — Xi + j- m —1 for all 1 ^ i ^ j ^ m — 1, i + j > m, 

mxi + i <F + m for all i G {1, . . . , m — 1}, 

mxk*+k* = F + m, 

Xi e Z for alH G {1, . . . , m — 1}, 

where k* — F (mod m). 

The following result follows from the above reasoning. 

Theorem 3.2. For fixed multiplicity m, counting the numerical semigroups for any Frobenius number F 
is doable in polynomial time. 

Remark 3.1 (Maximal Embedding dimension semigroups). A numerical semigroup is a maximal embed- 
ding dimension numerical semigroup (MED -semigroup for short) if its multiplicity equals its embedding 
dimension. Further details about these semigroups can be found in [T], |T7j . In [16] the authors give a similar 
characterization that the one used in the proof of Theorem \3.1\ but for MED -semigroups. Let MEDS(m) be 
the set of all MED -semigroups with multiplicity m G N\ {0}. Then, there is a one-to-one correspondence 
between MEDS(m) and the set of solutions of the following system of linear diophantine inequalities: 

Xi 1 for all i G {1, . . . , m — 1}, 

(4) Xi + Xj — Xi + j ^ 1 for all 1 ^ i ^ j ^ m — 1, i + j m — 1, 
Xi + Xj — Xi + j- m for all 1 ^ i ^ j ^ m — 1, i + j > m, 

Xi 6 Z for all i G {1, . . . , m — 1}. 

With this characterization and following the same scheme that the one used in the proof of Theorem \3.1\ 
we can state that for fixed multiplicity m, counting the MED-semigroups for any genus g is doable in 
polynomial time. Also, the analogous of Corollarv \3.1\ is true: Let g be a fixed genus, counting the number 
of MED-semigroups of genus g is doable in polynomial time. 
Tabled summarizes the results obtained up to genus 15. 
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Table 2. Number of MED-semigroups with given gender g and multiplicity m. 
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Analogous results may be stated for symmetric numerical semigroups or maximal embedding dimen- 
sion numerical semigroups, using the correspondence between the corresponding family and the polytope 
describing the Apery set. 

In the following section we exploit the use of short generating functions for the smallest multiplicities, 
where the polytopes are more manageable and formulas for the number of integer points inside them can 
be obtained for any genus or Frobenius number by using differente techniques. 

4. Numerical Semigroups with multipicities three and four 

In this section we analyze, using the above methodology, those numerical semigroups with multiplic- 
ities three and four. For these families we obtain the short generating functions defining the polytopes 
characterizing these semigroups with arbitrary genus or Frobenius number. Taking advantage of the 2- 
dimensional and 3-dimcnsional geometry of the corresponding bodies, we give formulas for the number of 
these numerical semigroups. 

4.1. Multiplicity Three. Let <S 3jS be the set of numerical semigroups with multiplicity m — 3 and genus 
g. The following result depicts the number of numerical scmigropus in S^_ g . It is based on the structure 
of the 2- dimensional polytope that characterizes the set iS 3jS . 

Theorem 4.1. The number of numerical semigroups of multiplicity 3 and genus g is — . 

Proof. Let P^, g be the polytope © when m = 3, that is, P^^ = {(x,y) : 2x — y>0, —x + 2y > —l,x + y = 
9,x>l,y> 1}. 

Since the integer points inside P^^ g must be in the line x + y = g we can reduce this polytope to the 
1-dimensional polytope P3 g = {x : | < x < }. Then, since y = g — x is integer if and only if x is 
integer, the number of integer points inside P 3 , s is #{P^ g n Z) = [^ip-J - [f ] = l 3 ^] ■ a 

The above result coincides with Corollary 10 in [TH] that was proved following a different strategy. 

Remark 4.1 (Maximal embedding dimension numerical semigroups with multiplicity 3). An analogous 
analysis can be done for MED-semigroups using the polytope characterization for these semigroups ([4|). In 

this case, the number of MED-semigroups of multiplicity 3 and any genus g > 2 is ^—^ — . 

Let £3,3 be the set of numerical semigroups of genus g and multiplicity 3 and S™ g the set of MED- 
semigroups of genus g and multiplicity 3. 

(1) The set iS 3;9 coincides with the set S™ g if and only if g = 2 (mod 3). 

(2) If g ^ 2 (mod 3), then, #iS™g = S™ g + 1, i.e., there exists only one numerical semigroup of genus 
g and multiplicity 3 that does not have maximal embedding dimension. 

Note that the above remark is not surprising. If S is an embedding dimension three numerical semigroup 
that is not of maximal embedding dimension (equals to 3 in this case) and genus g, then S = (3, n) for some 
n > 3. Since in this setting g(5) = F(s 2 )+1 = | (3 - l)(n - 1) = n - 1, we get that S = (3, g + 1). Eg = 2 
(mod 3), g + 1 is a multiple of 3, and S is not a numerical semigroup, and then there are no numerical 
semigroups with multiplicity 3 that are not of maximal embedding dimension. Otherwise (g = or 1 
(mod 3)) there is a unique numerical semigroup that is not of maximal embedding dimension. 

In the following we count numerical semigroups with multiplicity 3 and arbitrary Frobenius number. In 
this case, the set of numerical semigroups of multiplicity 3 and Frobenius number F is described by the 
system of diophantinc inequalities Q when m = 3. 
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Wc distinguish two cases, attending to the two possible choices for F: F = 1 (mod 3) and F = 2 
(mod 3). 

• If F = 1 (mod 3), then k* = 1, and the System ^ is reduced by substituting x\ by ^y^, the 



xi < and xi > £±1. Then, the number of possibilities for x\ are L"^ipJ — l""^"^! + 1 = 
L^-^-J — l~"^ir"] ^' coinciding with the case when F = 1 (mod 3). 

Proposition 4.1. T7ie number of numerical semigroups with multiplicity 3 and Frobenius number F (that 
is not a multiple of 3) is [^-^-\ — ["^p"] + 1- 

By using the same methodology, but applied to MED semigroups, we obtain the following result. 

Proposition 4.2. The number of MED semigroups with multiplicity 3 and Frobenius number F (that is 
not a multiple of 3) is L"^p"J — + 1- 

Remark 4.2 (Fixing genus and Frobenius number). With this procedure, adding to the polytope describing 
numerical semigroups with fixed genus, g, those constraints that fix the Frobenius number, F , we count the 
number of numerical semigroups with both values fixed. In the case when the multiplicty is 3, we obtain 
that all these polytopes ( varying F and g ) have a unigue integer point inside them. This fact reproves 
the well-known result that states that a numerical semigroup with multiplicity 3 is completely determined 
by its genus and its Frobenius number (see |18j). Actually, the unigue lattice point inside that polytope is 
(-£±2, 39 ~f~ 2 ) if F ee 1 (mod 3) or ( 3g ~f -1 , E ^-) if F = 2 (mod 3) being the unigue coordinates of the 
elements in the Apery set Ap(S,m) = {0,F + 3, 3g — F}. Then, the unigue numerical semigroup with 
multiplicity 3, genus g and Frobenius number F is S = (Ap(5*, m), m) = (3, F + 3, 3g — F). 

4.2. Multiplicity Four. Let <S4 iS be the set of numerical semigroups with multiplicity m = 4 and genus g. 
In the analysis for numerical semigroups with multiplicity 3, wc reduce the problem to dimension 1 since the 
polytope characterizing these semigroups is a 2-dimensional polyhedron and the equality constraint that 
fix the genus allows us to change one of the variables by a linear expression in the other. In this subsection 
we study, arguing analogously, the multiplicity four case. The main difference between this setting and the 
previous one is that in this case (m = 4) we are able to reduce the problem to a 2-dimcnsional polytope 
and generating functions tools are needed to count the integer points inside it. This is also the reason why 
to analyze larger multiplicities become much more complicated in 3 or more dimensions (see comment at 
the end of this section) . 

We distinguish different ranges of genus, since the geometry of the polytopes characterizing the numerical 
semigroups are slightly different. 



Proposition 4.3. Let P 4 , 9 = {(x, y, z) : 2x — y > 0, x+y — z > 0, — x+y + z > —1, —y+2z > —1, x+y + z = 
9,x>l,y>l,z>l}. 




• If g E [3,8] then #(P 4 , s nZ 3 



) = { 



4 



9 



3 



G 



7 



1 



if 9 
if 9 
if 9 
if 9 
if 9 
if 9 



3. 
4. 
5. 
G. 

7- 
8. 
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• If 9 > 9, then #(P 4 , 9 D Z 3 ) = #(T A (g) n Z 2 ) + #(T B (g) n Z 2 ) + #(R(g) n Z 2 ) - #(T c (.g) n Z 2 ) 

- T A (.g) = {(x,y) 6M 2 :3x + y>3,x< ^,y < §}, 

- T fl ( ff ) = {(x,y) £M 2 :x + 3y> fl -l,:c< < 

- R(g) = {(x,y) el 2 : <x< £±i,2^ <y < f}„ and 

- T c ( ff ) = {(x,y) eR 2 :x + y>3,x< ^.j/^ §}• 



Proof. The polytope P^g in R 3 can be projected in a lower dimension space by using the constraint 
x + y + z = g. By substituting each y by the expression g — x — z, we obtaing the following 2-dimcnsional 
polytope: P' = {(x, y) : 3x + y > g, z < |, x < , a; + 3y > 3 — 1, x + y < g — 1, x > 1, y > 1}. 

• Clearly, for 3 = 3, the number of integer points inside Pj^ is 1. 

• For g = 4, by using barvinok the generating function for the corresponding polytope is 

z\ z\z2 Z\Z2 



(l-f)(l-z 2 ) (l-|a)(l_ 2l ) (i_ Z2 )(i_ 2l ) 

By substituting in the above expression Z\ = £ and z 2 = £ 2 , and taking limit when £ — ► 1 we 
obtain that the number of integer points inside that polytope is 3. 
• For g = 5, barvinok output is 

Z 2 Z2 — zfz2 Z\Z 2 — ^2 Z 2 Z2 — ZiZ 3 , 



(1-%)(1- Zl ) ( i_fl )( i_ Z2 ) (l-z 2 )(l- Zl ) 

Proceeding as above we get that the number of integer points inside that polytope is 4. 
• For g = 6, invoking barvinok again, the generating function for the corresponding polytope is 

2 

— + Z 1 Z 2 —Z\zi,—z\z\ Z\z\ + z\z\ —Z2+Z1Z2 —Z\z\ — Z \Z2 



(l-f)(l-*i) (l-f)(l-z 2 ) (1-fXl-zx) (1 _f|)(i_^) (l-z 2 )(l-zi)' 

Repeating the above limit procedure, the number of integer points inside that polytope is 6. 
• For 3 = 7, the generating function is 

-Z 2 2 +Zl 3 Z 2 , ~Z 1 Z 2 3 - Z 1 4 Z 2 3 , Z! 2 Z 2 2 + Z! 3 Z 2 4 , -Zx 2 Z2 4 ~ Z 1 5 Z 2 



(1-^)(1- Zl ) (l_g)(l- 22 ) (l-^fl-Zl)' 

Being 7 the number of integer points inside that polytope. 
• For g = 8, barvinok yields 

-z 1 z 2 2 + z 1 4 z 2 ^~z 1 2 z 2 2 ~ zi 4 z 2 4 t z\z\ + z\z% i -ZlZ 2 5 + Zl 2 Z 2 2 ( -Zl 2 Z 2 5 - Zl 5 Z 2 



(i_a)(i_* 2 ) (i_g)(i_ Zl ) (1 _^ )(1 _ Z2 ) (i_ Z2 )(i_^)- 

And then, taking limit, as above, we obtain that the number of integer points inside that polytope 
is 9. 

• If 3 > 9, the polytope can be decomposed as the disjoint union P' n Z 2 = (Pa (3) H Z 2 ) U (Pb(3) n 
Z 2 ) U (R(g)\T c (g) n Z 2 ). Let (x,j/) G P' n Z 2 , if x < 2^1 j then, because of the constraint 
3a; + y > 3, (x,y) is in Pa (3) H Z 2 . If y < 2g ~ 3 , then, in view of x + 3y > 3 — 1, (x,y) is in 
Ts(g) nZ 2 . Otherwise, it is clear that the point (x, 3) is in R(g)\Tc(g). Note that since the vertex 
( , 2g ~ 3 ) is not an integer vector there are no integer points in any pairwise intersections of 
TA{g),Ts(g) and R(g)\Tc(g). Figure [1] shows this decomposition. 
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Figure 1. Decomposition of the polytopes for g > 9. 



□ 



In view of the above result, we can count the integer points inside each of the parts (triangles and 
rectangle) desctibing P^ g with g > 9, and adding all of them to obtain the number of points inside the 
whole polytope. By using the software barvinok for parametric polytopes (with parameter g > 9), we 
obtain the following formulas for the number of points inside each of the triangles. 

#(T A fo) n Z 2 ) = i - LSJ - 1) (29 - 3 Lf J - 2 Lf j - 2 - 3 L^J) , 



#(T B (g) n I?) = I (L^J - L^J) (-3 L^J + 2 LfJ - 3 L^J - 1) , 

#(Tc(g)<~)Z 2 ) = 1 (this is easily checkable analyzing the vertices of the rectangle and isosceles triangle Tc(g)). 
The generating function for a box H = [a, b] x [c, d) C M. 2 , is 



f(H;z 1 ,z 2 ) 



>1„M 



>1 J<*J 



(l- Zl )(l-z 2 ) (l-zi l )(l-z 2 ) (l- Zl )(l-z^) (l-zf 1 )(l-^ 1 )' 
Set Z\ = £ and z 2 = S, 2 , and take limit when £ — > 1. We obtain the well-known formula for the number of 
integer points inside the box H: 

(l+[b\ + \a])(l+[d\ + \b-]). 

By applying this formula to the rectangle R(g) we obtain the following descriptions for the numbers of 
points inside them. 

#(%)nz 2 



( 


5 + 1 




~2g + l~ 


♦0( 


9 




"2 5 -3" 


+0 




2 




8 


_2- 




8 



The following result summarizes the above reasonings. 



Theorem 4.2. Let S^ g be the set of numerical semigroups of multiplicity 4 and genus g > 9. Then the 
cardinality of S 4 . g is -g + § Lf J + +Lf J + \ J + L^J.9 ~ L^J Lf J - f L^J 2 - LfJ.9 + § LfJ 2 + 

LI J Lf J - 5 L^J + L^J Lf J - I L^J 2 + I L^J + I L^J 2 - L^J Lf J + L^J Lf J - L^J L Vl + 
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If we check the formula in the above result for those genus less than 9, we see that it is also verified for 
those genus g G {4, 5, 6, 7, 8}. For g = 3 is is known that there exists only one numerical semigroup with 
multiplicity 4 and genus 3. 



Corollary 4.1. Let S^ g be the set of numerical semigroups of multiplicity 4 and genus g. If g = 3, then 
S 4 , g = 1, otherwise the cardinality of S 4 . g is -g+ § [f J + + Lf J + H^irJ + L^J.? - L^J LfJ - |LhtJ 2 - 

Lf \9 + § Lf J 2 + Lf J Lf J - IL^J + L^J Lf J - f L^J 2 + IL^J + IL^J 2 - L^J Lf J + L^J Lf J - 

l^j r Vi + l^j - r Vi Lf j + l Vi r Vi - r Vi ■ 



In the following we describe the case when the multiplicity is four and the Frobenius number, F, varies. 
The set of numerical semigroups of multiplicity 4 and Frobenius number F is described by the of 
diophantinc inequalities in §3§ when m = 4. We distinguish three cases. 



If F = 1 (mod 4), then k* = 1, and the System ([3]) is reduced, after substituting x\ by -^-p: 



£2, £3 > 


1, 




4 x 3 - 4 x 2 


F + 


3. 


4x 2 +4x 3 > 


F - 


1. 


-x 2 + 2x 3 > 


-1 




4x 2 


<F- 


-2, 


4x 3 


< F- 


hi, 


x 2 , x 3 e Z. 







Using barvinok the number of solutions for the above parametric polytopc is 



F 2 -14F+141 
16 

-3 I F+l 2 
2 L 12 J 



1 I F+l 
4 L 12 _ 



-3 I F+l 
4 L 12 . 



F±5 J2 



BF 11 , F 
G 32 ' 32 



if 5 < F < 9, 
if 13 < F < 17, 
if F > 21. 



If F = 2 (mod 4), then fc* = 2, and the System Q is reduced, after substituting x 2 by ^-j^ 

xi,x 3 ^ 1, 
4x 3 -4x x ^ F + 2, 
4x 3 + 4x! > F-2, 
^-2^8x3^ F + 2, 

f + 2 8x1 SC2F + 6, 

xi, X3 e z. 

barvinok output for the number of solutions for the above parametric polytope is 





f + l 


4) 




-F + 5 


+ 




8 






8 
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If F = 3 (mod 4), then k* = 3, and the System ^ is reduced, after substituting X3 by ^-j^: 

xi,x 2 ^ 1, 
2xi — x 2 > 0, 
4xi +4x 2 > F + l, 
4x 2 - 4xi ^ F + 5, 

4xi <F + 3, 

4x 2 <F + 2, 

Xi, X2 6 Z. 

Invoking again barvinok, the number of solutions for the above parametric polytope is 



32 



7 

16' 



F 



19 
32 



F + l 


+ 


F + l 


2 1 


F 


F —- 

4 


F 


3 


F 


2 1 


F+l 


8 




8 


+ 4 


_12_ 


_12_ 


~ 2 


_12_ 


4 


8 



F. 



Theorem 4.3. TTie number of numerical semigroups with multiplicity 4 and Frobenious number F is 

. (^) 



1 F + 5 P 1 1 ^ + 5 
I "'"4 S 



F + l 



16 

5F 11 1 



F + 5 



+ ■ 



if F =1 (mod 4) and 5 < F < 9, 

if F = 1 (mod 4) and 13 < F < 17, 

if F = 1 (mod 4) and F > 21, 

i/F = 2 (mod 4), 



F^ , 7F _ 19 , 3 F+l 

32 t 16 



f[^J + [^J + i LSJ - 1 L-&J -ILSJ 2 -H^J F J/ ^ 3 (mod4) < 

otherwise. 



Remark 4.3 (Fixing genus and Frobenius number). It is well-known that fixing genus and Frobenius 
number is not enough to determine, in general, a numerical semigroup with multiplicity 4 (see Remark ^. 2\ 
for the analysis of the case of multiplicity 3). If we add to the polytope encoding the numerical semigroups 
with given genus, g, those constraints fixing the Frobenius number, F, we obtain this output in barvinok: 



F+3 
2 



4g-F-2 
4 



F+3 
2 



2g+F+5 
6 J 

2g+F+5 I 

6 J 

3F-4g+5 
4 

4g-F-7 
4 

2g-F 
2 

3F-4g+6 
4 
F-2 
4 

L^J 

3F-4g+7 

4 

2g+F+5 I 
6 I 

2g+F+5 
6 





ifF = 


1 


(mod 4) 


5F - 8.9 > 5, 


and 2g - F > 5, 


ifF = 


1 


(mod 4) 


Ag - F > 23, 


2g - F > 1 and2g-F < 3, 


ifF = 


1 


(mod 4) 


5F - 8.9 < 1, 


4.9 - 3F < 1 and 2g - F > 5, 


ifF = 


1 


(mod 4) 


4.9 - F > 11, 


4.9 - F < 19, 2.9 — F > land 2 5 


ifF = 


2 


(mod 4) 


8.9 - 5F < 2, 


F > 14, and 2g- F>2, 


ifF = 


2 


(mod 4) 


8.9 - 5F > 6, 


4.9 - 3F < 2 and 2g - F > 6, 


ifF = 


2 


(mod 4) 


8.9 - 5F > 6, 


F > 6, and2g-F < 4, 


ifF = 


2 


(mod 4) 


8.9 - 5F < 2, 


F < 10, and 2g - F > 2, 


ifF = 


3 


(mod 4) 


8.9 - 5F > 9, 


4g — 3F < 3 and 2g — F > 5, 


ifF = 


3 


(mod 4) 


8.9 - 5F < 5, 


and 2.9 — F > 5, 


ifF = 


3 


(mod 4) 


8.9 - 5F < 5, 


4.9 - F > 9, 2.9 - F > 1, and 2 5 


otherwise. 









Then, it is clear that different multiplicity four numerical semigroups can have the same Frobenius number 
and genus. If we analyze the three cases attending to the different possibilities for the embedding dimension 
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of the numerical semigroups we see that this result is not unusual. Let S be a numerical semigroup with 
multiplicity 4, genus g, Frobenius number F, and embedding dimension e(S). 

(1) If e(S) = 2, then S = (4 < 712}- From Selmer's formulas, 712 = and g — and then, S is 
completely determined by F, being S = (4, -^-j^)- 

(2) If e(S) = 3, then S = (4 < n% < TI3}. This case, is particulary difficult and is deeply studied in 
[19] . where the authors prove that S is completely characterized by F and g. 

(3) If c(S) = 4 (maximal embedding dimension case), then S = (4 < n2 < n$ < 714). From Selmer's 
formulas, 714 = F + 4 and g = n ^+ n ^+ F + 4 — ^Sp, being n 2 and n 3 all the integer solutions for 
r?2 + H3 = 4g — F + 2 with 4 < n 2 < n 3 < F + 4. This diophantine equation has, in general, 
more than one solution, actually, invoking barvinok for counting these solutions we obtain that 
the number of solutions is 

2 9 — L^-J — 4 l f 2g - F <2 and 4g - F > 9, 
2F-2g- [f J +2 if 2g — F > 3 and 4g - 3F < 3. 

However, in this system is not taken into account that gcd(n^, rij) = 1 for i,j = 1, . . . , 4, i ^ j, and 
then, the above formula is only a upper bound for the number of MFD -semigroups with multiplicity 
4, genus g and Frobenius number F . In case we want the exact number of these semigroups, we 
need to use the corresponding formula proceeding as in the general case, where barvinok output is 



F+l 
2 



2 



F+3 
2 



2g+F+l 

2g+F 
6 

3F+6 _ 
4 

2g+F+b 
6 

' gg+F+5 
6 



3F+7 
4 



if F 
if F 
ifF 

ifF 
if F 
ifF 
ifF 
ifF 



1 (mod 4 

1 (mod 4 

2 (mod 4 

2 (mod 4 

3 (mod 4 
3 (mod 4 
3 (mod 4 
5 and g = 



and -85 + 5F > 1 and 2g - F < 3 

and 8g - 5F > 3 and 4g - F > 15 and 4g - 3F < 1, 

and 8g - 5F < 6 and 2g- F >2 and 4g - F > 14, 

and 8g - 5F > 10 and 4.g - 3F < 2, 

and 8g - 5F < 1 and 2g — F > 5, 

and 8g - 5F < 1 and 4# - F > 13 and 1 < 2g - F < 3, 
and 83 - 5F > 5 and 4.g - 3F < 3 and 2g - F > 5, 
7, 



otherwise. 



The methodology applied in this paper for numerical semigroups with multiplicities 3 and 4 can be also 
applied for larger multiplicities. The case m = 5 has become much more tiresome, and the computations 
are available to the interested readers upon request. The quasy-polynomial for the number of numerical 
semigroups with multiplicity 5 and any genus has 9 pieces, being some of them particulary large. 
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